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Abstract. We study the asymptotics of sums of matricially free random variables, 
called random pseudomatrices, and we compare it with that of random matrices with 
block-identical variances. For objects of both types we find the limit joint distribu- 
tions of blocks and give their Hilbert space realizations, using operators called 'matri- 
cially free Gaussian operators'. In particular, if the variance matrices are symmetric, 
the asymptotics of symmetric blocks of random pseudomatrices agrees with that of 
symmetric random blocks. We also show that blocks of random pseudomatrices are 
'asymptotically matricially free' whereas the corresponding symmetric random blocks 
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are 'asymptotically symmetrically matricially free', where symmetric matricial free- 
ness is obtained from matricial freeness by an operation of symmetrization. Finally, 
we show that row blocks of square, lower-block-triangular and block-diagonal pseu- 
domatrices are asymptotically free, monotone independent and boolean independent, 
respectively. 



!>■ ! 1. Introduction and main results 

We have recently shown that the Hilbert space construction of the free product of 
states on C*-algebras given by Voiculescu [8] can be generalized to a framework which 
O ■ exhibits some matricial features [5]. 

We considered arrays of noncommutative probability spaces, for which we defined the 
matricially free product of states. The definition of this product is based on replacing the 
family of canonical unital ^representations of C*-algebras on a free product of Hilbert 
spaces by a diagonal-containing array of non-unital ^representations of C*-algebras on 
the matricially free product of Hilbert spaces. The crucial point is that products of 
these representations imitate products of matrices rather than free products, although 
the main features of the latter are still present. We studied the associated concepts of 
noncommutative independence called matricial freeness and related to it strong matri- 
cial freeness which can be viewed as a scalar-type generalization of freeness underlying 
other fundamental notions of noncommutative independence (monotone and boolean) 
and some of their generalizations (conditional freeness and conditional monotone inde- 
pendence) . 

More importantly, matricial freeness is closely related to random matrices. Their 
significance for free probability was discovered by Voiculescu [10], who showed that 
Gaussian random matrices with mutually independent entries are asymptotically free. 
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This connection was later generalized by Dykema [2] to non-Gaussian random matrices. 
Of course, the first indication that free probability might be related to random matrices 
was the central limit theorem for free random variables since its limit law was the 
semicirle law obtained in the classical work of Wigner [13] as the limit distribution of 
certain self-adjoint random matrices with independent entries. 

The main objects of interest in our study of matricial freeness and strong matricial 
freeness are diagonal-containing arrays of (in general, non-unital) subalgebras of a given 
unital algebra A, 

(A,j)(i,j)eJ with ACJC/x/ and A = : j E I}, 

equipped with an array of states (<fij) on A, with respect to which notions of noncom- 
mutative independence are defined. Here, the states ifiij and their kernels play a very 
similar role to that of one distinguished state ip and its kernel in free probability. Apart 
from square arrays, of particular interest are lower- (upper-) triangular arrays related 
to monotone (anti- monotone) independence introduced by Muraki [6]. In all arrays, the 
roles of diagonal and off-diagonal entries are quite different, which is a characteristic 
feature of this theory reminding the random matrix theory. 

Let (Xjj(n)) be an n-dimensional square array (or its diagonal-containing subarray) 
of self-adjoint random variables in a unital *-algebra A(n) which is matricially free with 
respect to the array (ifiij(n)) defined in terms of a family {(Pj{n))\<j< n . The sums 

S(n)= 

(i,j)eJ(n) 

where (J(n)) is an appropriate sequence of sets, called random pseudomatrices, remind 
random matrices, whereas convex linear combinations 

1 n 

/ v 

3=1 

where n G N, replace the states r(n) = E <g> tr(n), classical expectation tensored with 
normalized traces, under which distributions of random matrices are computed. 

Under suitable assumptions on the distributions of the variables (n) in the states 
(Pij(n), respectively, we can now study the asymptotic distributions of sums correspond- 
ing to blocks of random pseudomatrices in various states. These blocks are defined as 
sums of the form 

S P,q( n ) = Yl X ^ n ) 

(i,j)£N p xN q 

where [n] := {1,2, ... ,n} = Ni U N 2 U . . . U A^ r is a partition into disjoint subsets whose 
sizes grow proportionately to n (it is convenient to think of intervals). The dependence 
of the iVj's on n is supressed in the notation. In particular, we assume that the variances 
Vij(n) := ipi t j(n)(Xfj(n)) are identical within blocks. 

This leads to various asymptotic properties of sums of matricially free and strongly 
matricially free random variables. Let (pij(n) be defined according to 

<Pij(n) = <p(n) and <Pij(n) = (pj(n) for % ^ j, 
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where (p(n) is a distinguished state on A(n) and (<Pj(n))i<j< n is a family of additional 
states called conditions associated with (p(n) for any n G N. Then, as n — > oo, we 
obtain the following properties: 

(i) if the variances of square arrays are identical within blocks and rows, then the 
sums 

n 

A p (n) : = S Pjq (n), where 1 < p < r, 

9=1 

are asymptotically free with respect to (p(n), 

(ii) if the variances of block-lower-triangular arrays are identical within blocks and 
rows, then the sums 

p 

B p (n) := S Pt q(n), where 1 < p < r, 

9=1 

are asymptotically monotone independent with respect to <p(n), 

(iii) if the variances of block-diagonal arrays are identical within blocks, then the 
sums C p (n) := S PjP (n), where 1 < p < r, are asymptotically boolean indepen- 
dent with respect to (p(n). 

Let us add that non-asymptotic analogs of (i)— (ii) for finite sums of strongly matricially 
free random variables were proved in [5] and that (iii) also holds for finite sums of 
both matricially free and strongly matricially free random variables. Recall that the 
'strongly matricially free product of states', on which the definition of strong matricial 
freeness is based, is obtained from the matricially free product by restriction. However, 
it is worth to remark that the difference between blocks of matricially free and strongly 
matricially free random variables disappears asymptotically. 

The scheme of matricial freeness and its symmetrized version, in which ordered pairs 
are replaced by (non-ordered) sets consisting of one or two elements, called symmet- 
ric matricial freeness, is also used for arrays (ipij(n)) defined in terms of a family 
(<Pj(n))i<j< n of states on A(n) according to 

<Pij(n) = <pj(n) foranyi,j, 

and any n G N. Here, we do not a priori assume that the states (pj(n) are conditions as- 
sociated with a distinguished state. In the study of asymptotics, of importance become 
arrays defined by 'normalized partial traces' 

I ha 

q ieiv, 

according to ip Ptq (n) = ip q {n) for any p,q G [r] and n G N, where n q is the cardinality 
of N q . Then, as n — > oo, we obtain the following properties: 

(iv) if the variances are identical within blocks, then the array (S Ptg (n)) is asymp- 
totically matricially free with respect to (i^ Pjq (n)), 
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(v) if the variances are identical within blocks and form symmetric matrices, then 
the array (Z Ptq (n)) of symmetric blocks given by 

Z P,q( n ) : = Yl X i>M) 

where N pA = (N p x N q ) U (N q x N p ), is asymptotically symmetrically matricially 
free with respect to 

(vi) if the variances are identical within blocks of symmetric random matrices, then 
the array (T Piq (n)) of symmetric random blocks is asymptotically symmetrically 
matricially free with respect to (r p (3 (n)), 

where the array (t Pi9 (ti)) is defined by the family of partial traces r q = E <g> tr q (n) and 
tr q (n) denotes the normalized trace over the vectors indexed by N q . Here, by symmetric 
random blocks we understand symmetric blocks of symmetric random matrices in the 
approach of Voiculescu [10] and Dykema [2] . 

Concerning a non-asymptotic analog of (iv), we show in this paper that blocks 
(S Pj q(n)) are matricially free for any finite n under the stronger assumption that the 
variables have block-identical distributions. In turn, properties (v) and (vi) do not seem 
to have their non-asymptotic analogs. Nevertheless, they are the reason why we call 
the sums S(n) random pseudomatrices. On the other hand, let us also remark that 
we do not have an analog of (iv) for random matrices. Therefore, informally, random 
pseudomatrices can be viewed as objects which remind random matrices for large n if 
we consider blocks with symmetric variances, but exhibit different features if the vari- 
ances are not symmetric, which leads to triangular arrays and relations to monotone 
independence. 

In that connection we would like to mention the result of Shlyakhtenko [7] who estab- 
lished a connection between a class of symmetric random matrices called random band 
matrices and an operator-type generalization of freeness called freeness with amalga- 
mation. Our result (vi) gives a connection between another class of symmetric random 
matrices called symmetric random blocks with block-identical distributions and a scalar- 
type generalization of freeness called symmetric matricial freeness. At present, it is not 
clear to us whether there is a relation between these two approaches. 

In our previous work, we expressed the limit distributions of random pseudomatrices 
with respect to </?(n)'s and ^(n)'s in terms of certain functions on the class of colored 
non-crossing pair partitions, as well as in terms of some 'continued mult if r act ions' ]5]. 
Both these realizations showed that the limit laws can be viewed as matricial general- 
izations of the semicircle laws. Moreover, these distributions turned out to be related 
to s-free additive convolutions and the associated notion of freeness with subordination, 
or s- freeness [3], concepts motivated by the results of Voiculescu [11] and Biane [1] on 
analytic subordination in free probability. Using s-freeness, we also studied s-free and 
free multiplicative convolutions, which allowed us to establish relations between these 
convolutions (as well as some other multiplicative convolutions) and certain classes of 
walks on appropriately defined products of graphs [4]. 

In this paper, we give Hilbert-space realizations of the limit distributions of random 
pseudomatrices and their blocks under <p(n), ipk( n ) an d i>(n), where the underlying 
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Hilbert space in all considered cases is the matricially free product 

of the r- dimensional array (J~ Pt q) of Fock spaces, where the diagonal and off-diagonal 
Fock spaces are, respectively, free and boolean Fock spaces over one-dimensional Hilbert 
spaces. Using realizations on T and operators which play the role of Gaussian operators 
in our approach, as well as their 'symmetrizations', we prove (i)-(vi). 

In Section 2, we present the combinatorics of colored non-crossing pair-partitions. 
In Section 3 we recall the basic notions and facts concerning matricially free random 
variables. In Section 4 we introduce matricially free Gaussian operators on T . Using 
them, we find ^-realizations of the limit distributions for random pseudomatrices in 
Section 5. In Section 6 we prove that blocks of (strongly) matricially free arrays of 
random variables with block-identical distributions are matricially free with respect to 
a suitably defined array of states. In turn, asymptotic matricial freeness of blocks of 
random pseudomatrices is proved in Section 7, where we also find J-'-realizations of their 
limit joint distributions. In Section 8 we introduce the notion of 'symmetric matricial 
freeness'. In Section 9 we find the limit joint distributions of symmetric random blocks 
and their ^-realizations and we show that symmetric random blocks are asymptotically 
symmetrically matricially free. In Section 10 we obtain results on asymptotic freeness 
and asymptotic monotone independence of rows of pseudomatrices. 

2. COMBINATORICS 

The combinatorics of our model is based on the class of colored non- crossing pair 
partitions, to which we assign certain products of matrix elements whose indices depend 
on the colorings of their blocks. 

For a given non-crossing pair partition ir, we denote by B(ir), C(ir) and TZ(tt) the 
sets of its blocks, their left and right legs, respectively. If 7Tj = {l(i),r(i)} and ttj = 
{/(j),r(j)} are blocks of tt with left legs and and right legs r(i) and r(j), 
respectively, then Tii is inner with respect to Tij if < l(i) < r(i) < r(j). In that case 
7r j is outer with respect to 7Tj. It is the nearest outer block of iii if there is no block 
7Tfe = {l(k), r(k)} such that < l(k) < l(i) < r(i) < r(k) < r(j). Since the nearest 
outer block, if it exists, is unique, we can write in this case ttj = 0(71*), = o(l(i)) 
and r(j) = o(r(i)). If 7Tj does not have an outer block, it is called a covering block. It 
is convenient to extend each partition ix e A/"C^„ to the partition 7? obtained from ix by 
adding one block, say 7r = {0,m + 1}, called the imaginary block. 

Let F r (7r) be the set of all mappings / : B(ir) — > [r] called colorings of the blocks of n 
by the set [r] := {1,2,..., r}. Then the pair (tt, /) plays the role of a colored partition. 
Its blocks will be denoted 

B(7T,/) = {(7ri,/),(7r 2 , /),..., (TTfc,/)}, 

where we use pairs (tTj, /) since to each block 7Tj we shall assign entries of a matrix 
which depend on the colors of both 7Tj and 0(7^) for any i G [k]. If the imaginary block 
is used, it is convenient to assume that it is also colored by a number from the set [r]. 
In Fig. 1 we give an example with the notions defined above. 
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blocks 7ri={l,6},7r 2 ={2,3},7r 3 ={4,5} 
legs £(7r)={l,2,4}, 7?.(7r)={3,5,6} 
coloring /(tt 1 )=/, f{^ 2 )=i, /(7T 3 )=fc 
nearest outer blocks o(7r2)=7ri,o(7r3)=7T2 

FIGURE 1. A colored non-crossing partition 

Definition 2.1. Let (ir, f) be a colored non-crossing partition with blocks as above, 
where / G F r (n) and let B G M r (R) be given. For any < j < r we define 

where the functions fej : Z3(7r, /) — > R are given by the following rules: 

(1) &j(7Ti, /) = & Pi9 if f(m) = v and /(ofa)) = q, where < j < r, 

(2) bj(-Ki, f) = b p j if /(7Tj) = p and 7Tj does not have outer blocks, where 1 < j < r, 

(3) &o(7Tj, /) = b PtP if /(vrj) = p and 7Tj does not have outer blocks. 

The index j G [r] in &j(7r, /) can be interpreted as the color of the imaginary block 
(if j = 0, then the imaginary block is not needed). Finally, MC 2 m = for m odd and 
thus we shall understand in this case that the summation over ix G MC 2 m gives zero. 

When we sum these products over all possible colorings, we obtain numbers 

feF r (w) 

for any < j < r, which are used to express the limit laws of random pseudomatrices. 
Example 2.1. If we consider the partition it given by the diagram in Fig.l, we obtain 
&o(tt) = E 6 *AA« and b j( n ) = ^2hi b kA,j 

i,k,l i,k,l 

where j G [r] and all indices run over [r], which corresponds to all possible colorings 
from F r (ir) and is adequate for a square array B G M r (R). For other arrays, we obtain 
the same formulas except that range of the summation is smaller. 

Entries of the matrix B will be related to variances of matricially free random vari- 
ables obtained when computing their mixed moments. Therefore, it is natural to begin 
with an array (a it j) of random variables and assign a variable to each leg of the con- 
sidered partition. If the array is square, we assume that i,j G [r], whereas in other 
cases the pairs belong to a proper subset J C [r] x [r] which includes the diagonal. 
To a given non-crossing partition n G NC 2 m , where m is even, we can now associate a 
product of these variables in which indices can be interpreted as colors taken from the 
set [r] . These indices are related to each other in a natural way which refers to the way 
matricially free random variables can be multiplied to give a non-trivial contribution 
to the limit laws. The definition given below is based on this relation. 

Definition 2.2. We will say that the partition n G MC 2 m is adapted to a tuple of 
numbers (pi, q±, . . . ,p m , q m ) from the set [r] if 



imaginary block j 
I ; 1 



i k 

n n 

1 2 3 4 5 6 
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(1) (pi,qi) = (Pj,Qj) whenever {i,j} is a block of n, 

(2) qu = p (k) whenever k G TZ(n). 

The set of such partitions will be denoted A/"C^(pi, qi, ■ ■ ■ ,p m , q m )- 

If 7i G MC 2 m {p x , 9i, • • • ,Pm, qm), then the tuple (pi, qi, . . .,p m , q m ) defines a unique 
coloring of both n and n in which the block containing k G C{tx) is colored by pk for 
any k and the imaginary block is colored by q m . 

Example 2.2. Consider the partition ir shown in Fig.l and let pi,qi, . . . ,p 6 ,g 6 G [6] 
be given. Conditions (l)-(2) of Definition 2.2 lead to equations p\ = p§,P2 = Pz,Pa = Pb 
and ?2 — ?3 — Pi, 94 — 95 — Pi, 9i — 96) which gives four independent colors. Setting 
Pi — l,P2 — hPi — k and q$ = j, we obtain the coloring of Fig.l. Note in this context 
that Pi,P2,P4 may be interpreted as colors of the left legs of n and q Q as the color of 
the imaginary block. These equations lead to a product of variables of the form 

a pi,q6 a P2,Pl a P2,Pl a P4,Pl a P4,Pl a Pl m ) 

where each variable is taken from an array {a Ptq ) of matricially free random variables 
and is associated with one inner-outer pair of blocks of n. 

We close this section with the symmetrized version of Definition 2.2 which will be 
needed in a combinatorial formula for the limit joint distribution of symmetric random 
blocks. The symmetrization is obtained by replacing ordered pairs by subsets. 

Definition 2.3. We will say that the partition n G MC 2 m is adapted to a tuple of subsets 
({Pi, 9i}, • • • , {Pm, q m }) of the set [r] if 

(1) {pi, qi} = {pj, qj} whenever {i, j} is a block of n, 

(2) {p k , q k } H {p (k), q a (k)} ^ whenever k G K(ir). 

The set of such partitions will be denoted A/"C^({pi, qi}, . . . , {p m , 9m})- 

If 7r G NC 2 m {{Pi, 9i}, • • • , {Pm, qm}), the tuple ({pi, 9i}, • • • , {Pm, 9m}) defines color- 
ings of both 7i and n, called admissible, in which each block containing k G TZ(tt) is 
colored by p k or qk and the imaginary block is colored by p rn or q m . 

Example 2.3. Consider again the partition shown in Fig.l and let p±, qi, . . . ,pe, q& G [6] 
be given. Conditions (l)-(2) of Definition 2.3 take the form 

(1) {Pi,9i} = {P6,9e}, {^2,92} = {P3,9s} and {^4,94} = {P5,9s}, 

(2) {p 3 , q 3 } n {pi, qi} ^ and {p 5 , 9s} H {pi, qi} ^ 0. 

For instance, they are satisfied if (pi,qi) = (q e , Pe), (^2,92) = (93, P3), (P4,94) = (95, Ps), 
with qi = p i+ i for any i = 1, . . . , 5, which corresponds to the product of variables of the 
form 

a Pl ,P6 a P6 ,P3 a P3 ,P6 a P6 ,P5 Q P5 ,P6 °P6 ,Pl , 

where we have four independent indices Pi,P3,P5,P6, among which p 1 colors the imag- 
inary block and p 3 ,Pb,P% color the blocks of n (they are associated with the right legs 
of 7r). If we replace some indices from the set {pi,P3,P5,Pe} by the corresponding ^'s, 
we obtain other admissible colorings and the associated products of variables. 
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3. Matricially free random variables 

Let us recall the definition of matricially free random variables as well as the main 
results of [5], where we refer the reader for details. 

Let A be a unital algebra with an array {Aij) of not necessarily unital subalgebras 
of A and let ((fij) be a family of states on A. Here, by a state on A we understand a 
normalized linear functional. Further, we assume that each Aij has an internal unit ljj, 
for which it holds that aljj = l^a = a for any a G Aij, and that the unital subalgebra 
X of A generated by all internal units is commutative. If A is a unital *-algebra, then, 
in addition, we require that the considered functionals are positive, the subalgebras are 
*-subalgebras and the internal units are projections. 

However, the definitions given below are slightly more general than those in [5]. 
Namely, in contrast to the formulation given there, we do not assume any particular 
form of the considered array ((Pij)- Instead, we will distinguish the diagonal states as 
those of the form ipjj, where j G I, but we will not assume that they all coincide. This 
will enable us to use the concept of matricial freeness for a wider class of arrays, which 
turns out convenient in the formulation of our results. 

We shall use the subsets of (/ x I) m of the form 

A m = {((«i, l 2 ), (l2, h),---, («m, Wi)) : (ii,i2) 7^ («2, h) ^ ■■■ (im, Wl)}i 
where m G N, with their union denoted A = Um=i 

Definition 3.1. We say that (U,j) is a matricially free array of units associated with 
{Aij) and {(fiij) if for any diagonal state ip it holds that 

(1) ip{uiau 2 ) = (p{ui)ip{a)ip{u 2 ) for any a G A and U\,u 2 G X, 

(2) if a k G A ik j k n Ker<p ik j k , where 1 < k < m, then 

„ n , _ f ^{aa 2 . . .a n ) if ((i u j ± ), . . . , {i m , j m )) G A 
<p{al ilJl a 2 ...a m )-^ Q otherwise 

where a G A is arbitrary and {i\,ji) ^ ... ^ {i m ,j m )- 

Definition 3.2. We say that {Aij) is matricially free with respect to {(pij) if 

(1) for any a k G Ker(p ik j k (1 A ik j k , where k G [m\ and {h,ji) ^ . . . ^ {i m ,jm), and 
for any diagonal state ip, it holds that 

ip{a ± a 2 ...a m ) = 

(2) (ljj) is a matricially free array of units associated with {Aij) and {(pij)- 

Definition 3.3. The array of variables (ajj) in a unital algebra (*-algebra) A will be 
called matricially free {* -matricially free) with respect to {(p^j) if there exists an array 
of elements (projections) {lij) which is a matricially free array of units associated with 
A and {(fij) and such that the array of algebras (*-algebras) generated by a^j and ljj, 
respectively, is matricially free with respect to {(pij)- 

Slightly less general was the setting given in [5], where we assumed that all diagonal 
states coincide with some distinguished state ip and the off-diagonals states agree with 
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a family of additional states (<pj)j e i, called conditions associated with <p, which are 
defined by 

(Pj(a) = ^(cjabj) 

for some Cj, bj G Ajj H Kery? such that ip(cjbj) — 1 (if A is a *-algebra, Cj — b*). In this 
setting, it is natural to assume that (p and ^-'s are normalized according to 

<f(hj) = Sij and <Pj(l it k) = S j>k 

for any i,j, k. However, in this paper we will also use other arrays, for instance such in 
which all states in the j'-th column agree with some tfj, where (<Pj)jei is a given family 
of states on A which are not a priori assumed to be conditions associated with some 
distinguished state. This motivates the following definition. 

Definition 3.4. Let ip and (pj, j G /, be states on A and let ((fij), where G J and 
A C J C / x J, be an array of states on A. If ipjj = tp and ipij = ipj for any G J, 
we will say ((pij) is defined by ip and the family (ipj)j e i. In turn, if ipij = <pj for any 
G J, we will say that ((p^j) is defined by the family (ipj)j e i. 

In this context, let us remark that a Hilbert space setting, similar to that in [5], can 
be given for a family (<pj)j eI of product states on the unital *-algebra 

A-—U--A-- 

which are not defined as conditions associated with some distinguished state on A. 
Here, UijAij stands for the free product of C*-algebras without identification of units 
which is assumed to contain the empty word playing the role of the algebra unit. 

In this 'tracial' framework, if ('Hij,'Kij,£,ij) is the array of GNS triples associated 
with an array (Aij, <fij) of noncommutative probability spaces, the underlying product 
Hilbert spaces are of the form 

(oo 

where j & I and T-L^j is the orthocomplement of C^j in 1-Lij for each and their 
direct sum is 

^ = 0^J' 

jei 

with the canonical inner product used in all direct sums. The space % replaces the 
matricially free product of Hilbert spaces *ij(7~Lij,tiij) used in [5]. 

Nevertheless, H can be embedded in a larger matricially free product of Hilbert 
spaces, namely (H',g) = *fj~CH'ij,&j), where 



' 1 " < Hjj © if i — j anQ ^ if i = J 



and £jj is an additional unit vector for each j G I. In order to define appropriate 
product states on A associated with vectors £jj, we first trivially extend each cyclic 
representation Hjj : Ajj — > B(l-Ljj) to a non-cyclic representation 71^- ■ : Ajj — > B{H'j ■), 
and we keep the off-diagonal representations unchanged, namely 7r£- = 7Tjj for i ^ j. 
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Then, we take the free product A' = *fj^ij and observe that H is left invariant by A'(a) 
for any a G A. This allows us to define the unital ^representation 

X:A^B(H) as \ = \'\H 

which then leads to product states tpj on A defined by 

(fj = T)j o A, 

where rjj is the vector state on B(T-L) associated with where j G I. An equivalent 
formulation in terms of partial isometries leading to the product states tpj can also be 
given. 

The proposition given below provides the main motivation for allowing a more general 
class of arrays in the definition of matricial freeness. 

Proposition 3.1. The array of C* -algebras (Aij) viewed as *-subalgebras of A = 
UijAij is matricially free with respect to the array (ipij) of states on A defined by the 
family (<fj)j eI introduced above. 

Proof. In fact, it can be seen that an analog of [5, Proposition 2.3] holds for the states 
tpj := rjj o A defined above. However, in this case it suffices to take a k G Ai k j k , where 
k G [n] and ^ ... ^ (i n ,jn) (thus, we do not assume that (ii,ji) ^ ^ 

(i n ,j n )). If a k G Kenp iktjk for k G [n], then 

tp j (a 1 a 2 ...a n )=0 

for each j. Moreover, if a r = lj rJr and a m G Ker tp im j m for r < m < n, then 

i , _ ( <p j (a 1 ...a r _ 1 a r+1 ...a n ) if ((i r , j r ), . . . , (i n , j n )) G A 

<PMi-0-\ o otherwise 

Finally, for any a G A, u±, u 2 G X and k G /, it holds that 

ip j (u 1 au 2 ) = (p j (ui)(p j (a)(p j (u 2 ) and <Pj(li, k ) = S jtk . 

All these facts imply that (Aij) is matricially free with respect to the array (ipij), where 
A,j = <Pj for any ■ 

Remark 3.1. In this context, let us observe that if (Aij) is matricially free with respect 
to the array defined by tp and the associated conditions ((fj)jeij then (Aij) is not, in 
general, matricially free with respect to the array defined by (<Pj)jei since condition 1 
of Definition 3.2 does not need to hold if tp is replaced by tpj and we take ai,a m G Ajj. 

Let us assume now that for any natural n we have an n-dimensional square array 
of self-adjoint variables (Xij(n)) in a unital *-algebra A(n) equipped with an array 
of states (tpij(n)) defined by a family of states (<Pj(n))i<j< n - Blocks are defined by 
partitioning the set {1,2, ... ,n} into disjoint non-empty subsets, 

[n] = Ni U N 2 U . . . U N r 

where r G N, with dependence on n supressed in the notation, whose sizes increase as 
n — ?■ oo so that nj/n — > dj, where rij is the cardinality of Nj for any j G [r]. Then the 
numbers dj form a diagonal matrix D of trace one called the dimension matrix. 
Concerning the distributions of the considered arrays, we assume that 
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(Al) 
(A2) 



(X it j(n)) is matricially free with respect to (</?jj(n)) for any n G N, 
the variables have zero expectations, 



<Pij(n)(X itj (n)) = 



(A3) 



for all i,j G [n] and n G N, 

their variances are block-identical and are of order 1/n, namely 



u 



n 



(A4) 




for all i, j G [n] and n G N. 

In particular, if the distributions of the X^n^s in the states <Pi,j(n) are block-identical, 
assumptions (A3)-(A4) are satisfied. 

Following [5], where we studied limit distributions of random pseudomatrices under 
i/j(n), we consider now normalized partial traces 



where k G [r] and n G N. Easy modifications of the proofs given in [5, Lemma 6.1] and 
[5, Theorem 6.1] lead to combinatorial formulas for the limit distributions under partial 
traces given above. 

Theorem 3.1. [5] Under assumptions (Al)-(AJ,.), the limit distributions of random 
pseudomatrices under partial traces have the form 



where k G [r], m G N and B = DU, with D being the dimension matrix. Consequently, 
ip{n)(S m {n)) converges to J^weAfc 2 ^( 7r ) as n ~ * °°> where h{ji) = Y7k=i dkbk(n) . 

The above result refers to the 'tracial' framework which reminds the limit theorem 
for random matrices and is of main interest to us in this work. However, we have also 
shown in [5] that a similar result holds for the 'standard' framework which reminds the 
central limit theorem and involves the distributions of random pseudomatrices in the 
distinguished states (p(n). This can be phrased as follows. 

Theorem 3.2. [5] Under assumptions (Al)-(AJf.) for arrays of states ((pij(n)) defined 
by a distinguished state <p(n) and the associated conditions (v 9 j( n ))i<i<n for each n, 






where m G N and B = DU. 
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In the sequel we will derive Hilbert space realizations of the limit joint distributions 
of random pseudomatrices S(n) and their blocks under <p(n), ipk(n) and ifj(n). Interest- 
ingly enough, all these realizations are given on the same type of Hilbert space which 
is a matricially free product of Fock spaces. 

4. Matricially free Gaussian operators 

In this Section we shall introduce self-adjoint operators which play the role of matri- 
cially free Gaussian operators living in the matricially free product of Fock spaces. 

Recall that by the boolean and free Fock spaces over the Hilbert space H, respectively, 
we understand the direct sums 

oo 

F Q (H) = C£ © U and T(U) = C£ © H® m , 

m=l 

where £ is a unit vector, endowed with the canonical inner products. We shall use them 
to define an array of Fock spaces and their matricially free product. 

Definition 4.1. Let De an array of Hilbert spaces with distinguished unit 

vectors. By the matricially free product of we understand the pair (H,t;), 

where 



with U\j = Htj © C&j and f being a unit vector, with the canonical inner product. 
We denote it (H,0 *i'(H,. r C 



1,3) 



We already know that the matricially free Fock space, which is a matricially free 
analog of the free Fock space, is the matricially free product of an array of free Fock 
spaces [5]. Nevertheless, in order to find Hilbert space realizations of the limit laws 
of Theorems 3.1 and 3.2, it suffices to take a matricially free product of an array of 
Fock spaces, in which free Fock spaces are put only on the diagonal and the boolean 
Fock spaces elsewhere. Clearly, we obtain a truncation of the matricially free Fock 
space in this fashion. This structure is related to the difference between diagonal and 
off-diagonal random variables and is related to the difference between multiplication of 
diagonal and off-diagonal blocks of usual matrices. 



Definition 4.2. By the matricially free-boolean Fock space over the array % = (Hi,j) 
we shall understand the matricially free product 

ToCHij) if i^j 
and £jj denotes the distinguished unit vector in T^j. 



{T,i) = where T i4 



Remark 4.1. If we have a square array M. = ('Hij), where T-L^j = Hi for any i,j<El 
and (Hi) is a family of Hilbert spaces, then we have a natural isomorphism 
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since each tensor product 
is isomorphic to 

^r 1 ® ?c 2 • • • ® ^r™ 

for any «i 7^ i 2 7^ . . . 7^ i m and m, ni,ri2, ■ ■ ■ ,n m G N. Thus, J 7 is in this case also 
isomorphic to the strongly matricially free Fock space 71(H) introduced in [5]. 

Example 4.1. A simple example of a matricially free-boolean Fock space is the matri- 
cially free product of the two-dimensional array of the form 

(T \ _ ( h(Gi t i) l 2 (G lj2 ) \ 
[ ij) V HG2,i) k(G 2 , 2 ) J ' 

where Gjj = FS(gjj), the free semigroup on one generator gjj, where j G {1,2}, and 
Gij = Z 2 with the generator denoted g id for % 7^ j. Then J 7 = l 2 (G), where G is the 
'matricially free product of semigroups' Gij, by which we understand the subset of their 
free product *i,jGij given by the union 

00 

G = [j G in) 

n=0 

of disjoint subsets, where G^ consists of words of type g±g 2 . . . g n with g^ being 
an element of G® j , where Gr° ■ := Gij \ {e%,j}, where is the unit in G^, with 
((ii,ji), • • • , (i n ,jn)) G A and i n = j n . For instance, 

G (0) = {e}, 

G (1) = {g k 1A , g k 2 , 2 ■■ keN}, 

G {2) = {g 2 ,i9i,i, g li2 g 2 , 2 ■ k e N}, 

C (3) = {#2,202,101?1 ,9il9l,29™2 ,9l,292,l9™l ,02,101,2^2 : ^,^eN}, 

etc., with the remaining subsets consisting of words built from 'matricially free products' 
of powers of the generators, where the diagonal generators admit all natural powers, 
whereas the off-diagonal ones admit only the powers equal to one. 

Definition 4.3. Let A = (aiij) be a diagonal-containing array of positive real numbers 
and let (Hij) = (Ce it j) be the associated array of Hilbert spaces. By the matricially 
free creation operators associated with A we understand operators of the form 

Sij = ®i,jT*£(e id )T, 

where r : T — > ? r (@ i jHij) is the canonical embedding and the ^(e^'s denote the 
canonical free creation operators. By the matricially free annihilation operators and the 
matricially free Gaussian operators we understand their adjoints ^* ■ and sums denoted 
dj = %j + <hj, respectively. 

We shall assume now that A is a diagonal-containing subarray of a finite square array 
and that it is indexed by the set J (thus, A C J C I x /). Moreover, it is convenient 
to assume that A is the square root of another array B taken entrywise, i.e. 

a i,j = written A = \/~B 
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where B can be assumed to be a subarray of a square array. The dependence on A of 
the operators defined above is supressed in our notations. 

We shall prove below that the *-algebras Aij, each generated by and suitably 
defined unit ljj, respectively, are matricially free with respect to a suitably defined 
array of states on B^J 7 ). Namely, ljj is defined as the projection onto the subspace of 
T onto which the *-algebra generated by the creation operator acts non-trivially. 
To be more precise, let us introduce projections s it j and r it j for any j) G J which give 
an orthogonal decomposition of ly, i.e. ljj = r it j + Sjj and r^s^j = 1. Namely, Sjj is 
the canonical projection onto the subspace of J 7 spanned by tensors which begin with 
e it j for any G J. In turn, r it j is the canonical projection onto the subspace of T 
spanned by tensors which begin with for some k if i ^ j, whereas rjj is the canonical 
projection onto the subspace spanned by the vacuum vector and tensors which begin 
with e.,^ for k ^ j, where (j, k) G J. 

Two types of transformations on the considered operators will be performed: trunca- 
tions and symmetrizations. Here, we shall introduce truncated matricially free creation 
operators and truncated units by 

Pij = %3 P and Uj = UjP 

respectively, where P is the canonical projection onto T CQ. The truncated ma- 
tricially free annihilation operators and truncated matricially free Gaussian operators, 
respectively, will be denoted by p* ■ and u>ij = pij + p* ■ for any For uniformity, 
one can put P on both sides of the creation, annihilation and unit operators. 

Finally, in the case of finite dimensional arrays, which will be considered from now 
on, we use the same symbol to denote the sum of operators in a given array, like 

C= Yl 6j and u= Yl 

(i,j)eJ (i,j)eJ 

called the Gaussian pseudomatrix and the truncated Gaussian pseudomatrix, respec- 
tively. All adjoints are denoted in the usual way. Note that all these operators are 
bounded since all considered sums are finite. 

Proposition 4.1. Let ((pij) be the array of states on B^J 7 ) defined by the vacuum state 
if and the family where ipj is the vector state associated with ejj for any j G /. 

Then 

(1) the ip- distribution of Qj is the semicircle law of radius 2<x,j for any j, 

(2) the ipj -distribution of Q,j is the Bernoulli law concentrated at ±y/2ctij for any 
i 7^3, 

(3) the array (Cij) is matricially free with respect to (<Pij)- 

Proof. The first two claims follow easily from the definitions of the operators involved. 
Next, instead of proving the third claim, we will prove a slightly more general result 
that the array (Aij), where Aij = C(qj,^* j -, ljj) for any G J, is matricially free 
with respect to (ifij). Essentially, we proceed as in the free case [9], but we have slightly 
more complicated relations between the operators involved (see also Proposition 4.2 of 
[5]). In particular, one has to treat diagonal and off-diagonal subalgebras separately. 
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Let us consider first the off-diagonal case, when the creation and annihilation operators 
satisfy relations 

stj%j = b i,3 r i,j and ^ = 0, ^ = 
for any % ^ j. In that case we have additional relations 

r i,j^i,j = 0> ^i,j r i,j = S i,jSi,j — ^i,j s i,j = 0- 

Using all these relations and their adjoints, as well as equations 

^(nj) = l, i/>j(sij) = 

we deduce that an arbitrary noncommutative polynomial from Aij H Ker(^), where 
i 7^ j, is spanned by qj, and Sjj. In the diagonal case, the situation is similar 

to that in the free case since each pair of diagonal creation and annihilation operators 
satisfies the relation 

= h j-j l j-j with <P(hj) = 1, 

and therefore any polynomial from Ajj H Ker(</?) is spanned by ljj and ^fyj, where 
p + q > 0, for any j. Moreover, s^q^i = and q k ,iSij = ^k,i for any ^ (k,l). 

Hence, to prove matricial freeness of the array (Aij) with respect to (<fij), it suffices 
to show that 

for suitable powers Pi,qi, ■ ■ ■ ,Pm, a m that depend on whether the corresponding ope- 
rators are diagonal or not, since tp(si 1 j 1 . . . s« m j m ) = for any off-diagonal pairs ^ 
. . . 7^ (i m ,jm)- At this point we can use the same inductive argument as in the free 
case [9], which gives the above 'freeness condition'. Moreover, the definition of each 
shows that it is the projection onto the subspace onto which the *-algebra generated 
by acts non-trivially, which implies that the array (ljj) is the matricially free array 
of units. ■ 

Note that in Proposition 4.1 each ipj is a condition associated with ip since there exists 
bj G Ajj H Kenp, namely bj = Sjj/oejj, such that ipj(a) = <p(b*abj) for any a G B(F). 
However, a similar result is obtained for arrays of states defined by the family (ijjj)j e i, 
where ipj is now the vector state on BiT Cfl) associated with ejj for any j G I. 

Proposition 4.2. Let (ipij) be the array of states on B^FQCQ) defined by the family 
(ipj)j € i, where ipj is the vector state associated with ejj for any j . Then 

(1) the ipj -distribution of ujjj is the semicircle law of radius 2ajj for any j , 

(2) the if) j- distribution of Uij is the Bernoulli law concentrated at ±y/2otij for any 

(3) the array (wjj) is matricially free with respect to ijpij). 



Proof. The proof is similar to that of Proposition 4.1. 
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5. FOCK-SPACE REALIZATIONS OF LIMIT DISTRIBUTIONS 

Using the matricially free Gaussian operators and their truncations, we will find 
realizations on T of the limit distributions of Theorems 3.1-3.2. 

To each tt G MC 2 m for m even we assign natural products of creation and annihilation 
operators. First, to each tt G NC? m we assign the sequence e(7r) = (ei, €2, ■ ■ ■ , e m ), where 
€j = 1 whenever j G 7?.(7r) and e& = * whenever k G £(7r). Then, to each tt G A/"C^ we 
assign products of creation and annihilation operators 

?(tt) = ? £ V 2 . . . <f m and p^ir) = p ei p €2 ...p 6m 

for any j G [r], where (ei, e 2 , . . . , e m ) = e(vr). 

There is a nice relation between the expectations of these products and numbers bj (tt) 
defined in Section 2. In the case of b (tt) we will use the expectations in the vacuum 
state ip and for the remaining j's we take ^-'s associated with vectors ejj. Finally, to 
obtain b{j), we will use the convex linear combination of states of the form 

r 

i=i 

where numbers di, d,2, ■ ■ ■ , d r are taken from the dimension matrix D. 

We are ready to give Fock-space realizations of the limit distributions of random 
pseudomatrices in terms of the distributions of Gaussian and truncated Gaussian pseu- 
domatrices. 

Lemma 5.1. For the limit distributions of Theorem 3.2, it holds that 

7T€AfC 2 m 

for any m G N, where ( is the Gaussian pseudomatrix. 

Proof. If m is odd, both sides of the first equation are clearly zero. Therefore, suppose 
that m is even. We have 

y?(C)= E ¥>(^V a ...?*")= £ 

ei,e2,...,e m g{l,*} 7reA/"C^ 

where we used the fact that products of creation and annihilation operators corre- 
sponding to sequences (ei, 62, . . . , e m ) which are not associated with non-crossing pair 
partitions give zero contribution (this follows from the definition of the ■). The as- 
sertion of the lemma follows then from the combinatorial formula 

for each it G MC 2 m . To prove it, first observe that 

f^Pugi ■ ■ ■ 4Z,gJ = unless (ei, • • • , e m ) = e(ir) 

for some tt G AfC^. Now, if tt G AfC 2 m (pi, qi, ...,p m , q m ) for some p u q u . . . ,p m , q m G [r] 
and (e 1? . . . , e m ) = e(7r), then we claim that 
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where / is the coloring of n defined by p k , k G C{tc). This claim can be proved 
by induction. Suppose that the last annihilation operator in this mixed moment is 
indexed by k, i.e. e k = *. Then {k, k + 1} must be a block which does not have any 
inner blocks. The corresponding product of operators is of the form ^i k j k ^i k j k - Two 
cases are possible: 

(1) If k = m - 1 and p k = q k , then ^ k:q ^ Pkm acts on Q and gives b Pk ^ k Vt. 

(2) If k < m — 1 and q k+1 = p k +2, then ^ km ^ Pk<qk acts on some simple tensor h and 
gives b Pkt(lk h. Here, q k colors the nearest outer block of {k, k + 1}. 

If we repeat this procedure for the product of operators corresponding to the partition 
7r' obtained from n by removing block {k, k + 1}, we obtain the product of b p ^s which 
appears in the combinatorial formula for bo(7T,f) after a finite number of steps. If 
we fix 7r and sum over all p±, qi, ■ ■ ■ ,p m , q m to which ir is adapted (Definition 2.2), we 
obtain in fact the summation over p k ,k G C(tt), equivalent to the summation over 
F r (n), which proves the formula for bo(n) since mixed moments corresponding to those 
Pi, q±, . . . ,p m , q m to which n is not adapted are equal to zero. This completes the proof. 
■ 

Lemma 5.2. For the limit distributions of Theorem 3.1, it holds that 

for k G [r] and m G N and hence J^-jreAfc 2 K 71 ) = ^( wm ); where u> is the truncated 
Gaussian pseudomatrix. 

Proof. The proof is similar to that of Lemma 5.1 and is based on the analogous 
combinatorial formula 

b k (ir) = ip k (<p(it)) 

for any 1 < k < r and 7r G MC 2 m) where m is even and positive. The proof of that for- 
mula reduces to showing that if (ei, . . . , e m ) = e(ir) for some n G J\fC m (pi, qi, . . . , p m , q m ) 
and k = q m , then 

Mp e P \, qi ---p e P 2, q J = h(ir,f) 

where / is the coloring of n defined by p k , k G C(tt), with k coloring the imaginary 
block. As compared with the proof for ip, instead of acting on Q, we need to act on 
e k k . However, thanks to the projection P onto JGCO, which appears in the definition 
of the p e k ^ k , each e k k plays the role of the vacuum vector with respect to the p\ k for 
any i, k, including the case when i = k, since p ktk e k k = 0. Therefore, the arguments 
are similar to those for </?, except that pi )k e k)k is non-zero for arbitrary i. In terms 
of diagrams, this means that each covering block of n contributes bi tk for various i's. 
Finally, when we use the definition of ip and sum over k G [r], we obtain for each k 
the extra factor d k which appears in the combinatorial formula for b{it). Consequently, 
b(7r) = ip(p(ir)) for any such re, which then leads to the formula for ^(u n ). ■ 

Example 5.1. Consider the moment corresponding to the partition ir G MC\ consisting 
of two blocks: {1,4} and {2,3}. Let A = \fB G M 2 (R) be the square root taken 
entrywise, where B = DU and we set ai 5 i = a, = (3, «2,i = 7, «2,2 = 5- Then 

b (n) = a 4 + a 2 -f 2 + 5 2 (3 2 + 5 4 . 
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On the other hand, 

q 2 Q = a 2 (ei 5 i ® e^i) + a7(e 2 ,i ® e^i) 
+ /35(ei, 2 <E> e 2>2 ) + <5 2 (e 2)2 ® e 2 , 2 ) 
and thus the 'symmetric' action of (?*) 2 gives exactly bo(n)Q. 
Example 5.2. For the same partition n as in the above example we obtain 



A similar expression is obtained for d 2 ip2{p* P* pp) with a interchanged with S and (3 
interchanged with 7. The sum of both expressions agrees with b{ji). 



We show in this section that blocks of finite-dimensional arrays of matricially free 
random variables are matricially free with respect to an appropriately defined array 
of states. This is the analog of the property of free random variables which says that 
families of sums of free random variables are free. 

Definition 6.1. Suppose that the n-dimensional array pQj) of variables from a unital 
algebra A is matricially free with respect to {<fij) and let [n] — Ni U 7V 2 U . . . U N r be 
a partition of [n] ino disjoint non-empty subsets. The sums of the form 



where p,q G [r], will be called blocks of the pseudomatrix S = Ylij-^-ij- We will 
say that the array (X it j) has block-identical distributions with respect to ((pij) if the 
(/^-distributions of X it j are the same for all G N p x N q and fixed p, q G [r]. 

We need to define the associated array of block units which satisfy the conditions of 
Definition 3.1. In order to construct them on the level of noncommutative probability 
spaces, let us first return to the more intuitive framework of the matricially free product 
of representations of the array (Aij) of unital C*-algebras on the matricially free product 
of Hilbert spaces 



= rf 1 (a 4 + aV + /3V + 7V) 



On the other hand, 

V>(p(tt)) = di^i(p(7r)) + d 2 Mp(^)) 
where p(n) = p*p*pp, and we get 

P 2 ei,i = 7 (<K e 2,2 ® £2,1 <S> ei,i) + (3(e lt2 <S> e 2 ,i <8> e X)1 )) 
+ a (a(ei,i ® ei,i ® e^i) + 7(e 2 ,i ® ei,i ® e X)1 )) 

which, in view of the 'symmetric' action of the adjoint, gives 

d,MP*P*PP) = di(l 2 5 2 + /?V + a 2 7 2 + a 4 ). 



6. Matricial freeness of blocks 




(i,j)£N p xN q 



(n,0 
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and then carry over the corresponding definition to the algebraic framework of noncom- 
mutative probability spaces. We assume that each Aij is equipped with an internal 
unit ljj and a state ipij. It is the matricially free product of Hilbert spaces on which 
one defines the canonical ^representations of (Aij,tpij). Namely, let (Hij , TTi,j , be 
the associated GNS triples, so that <fi,j(a) = (^i,j(o)Ci,jjCi,j) f° r an y a e By \j 
we denote the canonical ^representation of (Aij,(pij) on (%,£). Using these represen- 
tations and appropriate partial isometries, we have defined in [5] their matricially free 
product A = *fjTTi,j which maps UijAij into B(T-L). 

In this connection note that A does not send the units ljj onto the unit of B{T-L). In 
fact, A(ljj) = r itj + Sjj, where 

rij = projection onto T-L(i,j) 
Sij = projection onto K,(i,j) 

where K,(i,j) = H^j <S>T-L(i,j) and T-L(i,j) denotes the subspace of T-L onto which the 
left free action of X(Aij) is non-trivial. By the left free action of X(Aij) we understand 
the left action onto the subspace spanned by simple tensors which do not begin with 
vectors from T-Q^ and, in addition, by £ if i = j. Clearly, _L s it j for any fixed i,j and 
thus their sum is the canonical projection onto the subspace of H onto which X(Aij) 
acts non-trivially. 

In the propositions given below we state basic properties of the projections involved 
that are needed for the construction of block units. Note that the case of card(J) = 1 
is trivial from the point of view of matricially free product structures and that is why 
it is not treated. By C[aj, i G I] we denote the unital algebra of polynomials in the 
commuting indeterminates a^i G /. Recall that X stands for the commutative unital 
algebra generated by the units A(l ij ) = 

Proposition 6.1. If card(J) > 1, then X = C[r it j, Sij,p^ : i,j G I], where p^ is the 
canonical projection onto C£. 

Proof. Observe that we have relations 

j,j = = s j,k and ljjlfcj = Sjj 

whenever i ^ j ^ k, which implies that Sj^,p^ G C[ljj : i,j G I] for any j, k, and since 
rj,k = h,k - s j,k, it holds that C[rij,Sij,p^ : %, j G I] C C[l id : i,j G /]. The reverse 
inclusion is obviously true, hence the proof is completed. ■ 

Proposition 6.2. If card(J) > 1 and Ji, J 2 C I are identical or disjoint, then the 
algebra C[ljj : i G J\, j G J 2 ] contains the canonical projection lj lt j 2 onto the subspace 
of H onto which the algebra generated by {X(a it j) : a it j G Aij, i G J\,j G J 2 } acts 
non-trivially. 

Proof. We want to construct a projection lj^jj onto the subspace of H that would 
be suitable for the left action of the A(djj), where % G Ji and j G J 2 . If Ji fl J 2 = 0, 
then this subspace is the orthogonal direct sum 

«Ji,j a = ^,i)©0(/c(i,i)e^(j,j)eco 
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and the associated canonical projection is 

1 Ji,J 2 = J2 S id + ^( S hi + r 3d-n) 
ieJi,jeJ2 ieJ2 

ieJi,jeJ2 

where k G J\ is arbitrary. In turn, if J\ = J 2 = J, then we obtain the orthogonal direct 
sum 

= C£ © (/c(j, j) © j) e co 

and the corresponding canonical projection takes the form 

jeJ 

where indices A;, / are such that k ^ I and otherwise are arbitrary elements of J. We 
have used the fact that = l^lj-j for any i ^ j. This completes the proof. ■ 

The proof of the above proposition enables us to construct block units which are 
internal units in the commutative algebras C[S p ^ q , l p , q ], where p, q G [r], each generated 
by a block and the corresponding block unit. Namely, we set 

lp, ? := lN p ,N q 

for any p,q G [r], where the right hand side is defined in terms of ljj's in exactly the 
same way as in the proof of Proposition 6.2. 

It will also be useful to introduce the following terminology. Namely, if 

a k e A ikJk n Ker(ip ikJk ) for 1 < k < n and ((i 

1 ) Jl J ) • • • i l^nj in)) e A, 

we will say that the product a x a 2 ... a n is in the matricially free kernel form with respect 
to (ip i:j ). 

Theorem 6.1. Let ((pij) be the array of states on A defined by the family (<Pj)i<j<n> an d 
let (ip Pt q) be defined by the family of associated normalized partial traces ijp q )i< q < r - If 
(X it j) is matricially free (strongly matricially free) and has block-identical distributions 
with respect to (</?ij), then (S p>g ) is matricially free with respect to (ip Pt q)- 

Proof. We will assume that (X it j) is matricially free since the proof for strong matri- 
cial freeness is analogous. Clearly, the unital algebra generated by (l p , g ) is commutative. 
We claim that the array (l p , ? ) is a matricially free array of units associated with (S Ptq ) 
and (i[)p t q). The proof of condition (1) of Definition 3.1 for (ip Ptq ) follows easily from the 
same condition for ((pij). To prove condition (2), we need to evaluate mixed momeints 
of type 

<Pj(wl p>q WiW2...W m ) 

where G Kei(ip Pkjqk ) is a polynomial in S PkSk for each k G [to] and j G [n] and the 
product W1W2 ■ ■ -w m is in the matricially free kernel form with respect to (ip p , q )- We 
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shall reduce the computations to moments of type ipj(wl i; jViV 2 ■ ■ ■ Vh), where v±v 2 ■ ■ ■ Vh 
is in the (strongly) matricially free kernel form with respect to ((pij)- For that purpose, 
let us express each power of S PkAk which appears in in terms of variables which are 
in the kernels of the tpij. This procedure, described in more detail below, is applied to 
w m , w m _i . . . , Wi (in that order). 

Let w(X) = J2 s r=0 c r X r be an arbitrary polynomial and let S p>q = J2^j) eNpXNq Xij. 
We decompose each positive power of which appears in w{S p ^ q ) as 

= Wj) + ^«Wj)i«. 

where i G N p ,j G N q . Using the fact that is a matricially free array of units, we 

observe that in the computations of the mixed moments of the given type we can use, 
without loss of generality, polynomials of the form 

s 

w {Sp,q) = Yj Yj Yj C h, , ---,irW+i^h,il) • • • C^l,ir+i) + C P>1 Yj 

r=l h,...,i r+1 ni,...,n r (i,j)eN p xN q 

with summations over ii, . . . , i r+ \ and ni, . . . , n r run over some finite sets of natural 
numbers, with n\ + . . . + n r < deg(u>), since the remaining terms will give zero con- 
tribution to the considered moment if the product of variables standing to the right of 
this polynomial is in the matricially free kernel form with respect to (<Pi,j), which is the 
case if we carry out our computations going from the right to the left. 

Note that if p ^ q, then s — 1, but if p = q, then s < deg(w). Let us also point out 
that thanks to our assumption that the array {X it j) has block-identical distributions 
with respect to (<fij), the same constant c P)(? stands by each ljj for any i G N p ,j G N q , 
i.e. it does not depend on Moreover, 

Vw( l ^ = — ^ U ^ = ^Y^ 1 ^ = n P 

(i,j)£N p xN q q k€N q (i,j)£N p xN q q i,j 

for any p, q, whereas the first sum in the above expression for w(S Piq ) belongs to 
Ker('0 Pi9 ). These arguments lead us to the conclusion that in the computations of 
mixed moments of the considered type we can take each polynomial Wk to be of the form 
w(S Pktqk ) given above, with c Pktqk = for each k G [m] since the product W1W2 ■ ■ .w m is 
assumed to be in the matricially free kernel form with respect to (<fii,j)- 

Therefore, each moment of type ipj(wl M wiW2 ■ ■ ■ w m ) is a sum of moments of type 
(pj(wl P)q v it>2 • • • Vh), where the product v\v 2 ... is in the matricially free kernel form 
with respect to Similarly, each moment of type ipj (wwiw 2 . . . w rn ) is a corre- 

sponding sum of moments of type ipj{wv\v 2 . . .Vh) since the reduction described above 
does not depend on what stands before the product wiw 2 . . . w m . It remains to observe 
that under each tpj, the block unit acts as the projection onto the linear span of 
products ViV 2 . . .Vj which are in the matricially free kernel form with respect to (tpij) 
and begin with V\ G A^^, where i\ G N q . The proof of that fact follows from the 
definition of 1 PA expressing it in terms of units 1^ which, under ipj, act as projections 
onto the linear span of v±v 2 . . . v m which are in the matricially free kernel form with 
respect to (<fij) and begin with v\ G where %\ = k. Of course, we use here the 

fact that (p,q) ^ This completes the proof of condition (2) of Definition 3.1. 
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Moreover, using the normalization conditions for ljj's, we obtain the normalization 
conditions 



which completes the proof that is a matricially free array of units. 

Finally, the proof of condition (1) of Definition 3.2 is similar to that of condition 
(2) of Definition 3.1 presented above and is based on reducing the computations of the 
mixed moments ipj(wiW2 ■ ■ ■ w n ), where W1W2 ... w„ is in the matricially free kernel form 
with respect to {ip p , q ), to mixed moments of type tpj{v\V2 ■ ■ ■ Vh), where v if 2 ... is in 
the matricially free kernel form with respect to ((fiij)- This completes the proof. ■ 



In this section we study the asymptotic joint distributions of blocks of random pseu- 
domatrices. In particular, we show that they are 'asymptotically matricially free', which 
is a notion analogous to asymptotic freeness. This generalizes the results of Section 6, 
where we proved matricial freeness of blocks in the case when the array of matricially 
free variables has block-identical distributions. 

For that purpose we will use a realization on the Fock space J 7 . Having established the 
realization of the limit laws of random pseudomatrices S(n) under normalized partial 
traces in Lemma 5.1, it is natural to expect that it can be carried over to the level of 
blocks of S(n) of the form 



for any n, where we require that the arrays (Xij(n)) of self-adjoint random variables 
satisfy the assumptions of Theorem 3.1. 

Before we proceed with examining the limit joint distribution of these blocks, we 
define the notion of asymptotic matricial freeness, following the analogous notion of 
asymptotic freeness. By the block units we shall understand units l p>q (n) defined for 
each n in terms of ljj(n)'s in exactly the same way as in Section 6. In addition to the 
states f(n) and ip(n) on the algebras A(n), we will also use normalized partial traces 
ip q {n), where q G [r] and n G N, defined as before in the case of fixed n. Informally, 
the states ip q (n) will play the role of conditions which converge to the conditions ip q of 
Section 4 as n — > 00. 

Definition 7.1. Let {i] Ptq {n)) be an r-dimensional array of functionals on the algebra of 
noncommutative polynomials C{S Ptq (n), l Pjq (n),p, q G [r]). We will say that (S Piq (n)) is 
asymptotically matricially free with respect to (r) p>q (n)) if these functionals have point- 
wise limits {r] p , q ) as n — > 00 with respect to which the limit array is matricially free. 

Theorem 7.1. Under the assumptions of Theorem 3.1, the joint ip q (n) -distributions of 
blocks and block units converge to the joint ip q - distributions of the truncated matricially 
free Gaussian operators and the truncated units, respectively, as n — >■ 00. 

Proof. This result is a refinement of Theorem 3.1. The combinatorial arguments 
referring to non-crossing partitions used in [5, Lemma 6.1] can be repeated except that 




7. Asymptotic matricial freeness of blocks 



(i,j)£N p xN q 



ASYMPTOTIC PROPERTIES OF RANDOM MATRICES AND PSEUDOMATRICES 



23 



since we now take products of the S Ptq (nys instead of powers of S(n), all indices run 
only over some of the subsets N\,N2,...,N r of the set [n]. Moreover, it suffices to 
consider the case of m even, say m = 2s, since if m is odd, both sides are zero, by 
standard arguments. We will first show that 

lim il> q (n)(S pu<n (n) . . . S Pm , Qm {n)) = ip q ( 

n—too 

for any q ^ 0. Thus, for given q and p±,qi, . . . , p m , q m , in order to compute the left-hand 
side of the above equation, it suffices to take into account only the mixed moments 

vM){X ildl (n) . . .X imdm (n)) 

in which (i k ,j k ) G N Pk x N qk for any k G [m\ and j G N q . 

However, in the limit n — >■ oo, further reductions takes place. As in the case of ip(n)- 
and y?(n)-distributions of random pseudomatrices investigated in [5, Lemma 6.1], the 
sum of mixed moments of matricially free random variables in the states <fj(n) which 
correspond to a partition which is not a non-crossing pair partition is 0(1/ 'yjn). 

Moreover, it suffices to take into account the mixed moments of the above type 
which are associated with n G NC 2 m {j>\, qi, ■ ■ ■ ,p m , Qm) for q = q m since the remaining 
moments vanish. However, these moments are the same for all ■ ■ ■ ,i m ,j m ), for 

which 7r G MC 2 m [i\,2\i ■ ■ ■ , im,jm) and j m = j since the variances of (Xij(n)) are block- 
identical, namely fjj(n) = u Ptq /n whenever G N p x N q . Thus we can use the 
matrix elements of U — (u Ptq ) to compute our mixed moments and express them as 

^(n)(X, lun (n) . . .X Wm (n)) = 

where / is the unique coloring of n G MC 2 m (j)\, qi, . . . ,p m , q m ) defined by indices p k , 
where k G C(tt), with q coloring the imaginary block. Since the coloring / is uniquely 
determined by p 1 , q 1 , . . . ,p m , q m , we will use the simplified notation u*(ir) = u q (n, f). 

For such Ti it then remains to enumerate the tuples (ii,ji,---,i m ,j m ), for which 
7r G MC 2 m [i\,2\i ■ ■ ■ i im, jm)- Since the coloring of n defined by these tuples is determined 
by independent indices associated with Pi(i), ■ ■ ■ ,Pi( s ), where 

£(ir) = {i(l),/(2),... ,/(«)}, 

and by q = q m , this enumeration boils down to enumerating the indices . . . ,i^ s ) 
and j = j m . An appropriate inductive argument for this fact can be easily provided. 
Instead of giving a formal proof, we refer the reader to Example 2.2 and to a similar 
argument given in the context of random matrices (see the proof of Theorem 9.1). 
Therefore, for given n, the cardinality of the corresponding set of equal mixed moments 
is of the same order as 

©n(vr) = n qm n Pl(1) . . . n Pl(s) = 0(n s+1 ) 

as n — > oo, where q m = q. In fact, the exact cardinality may be slightly smaller 
than Q n (n) since all independent indices i[, where / G C(tt), must be different in this 
computation to give a non-zero contribution to the limit. In the formula for n (7r) given 
above this is not the case since the mapping C(tt) — > [r] given by k — > p(k) may not be 
injective. Nevertheless, we can substitute O n (7r) for the cardinality of the considered 
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set when taking the limit n — > oo since the difference between these two cardinalities 
is 0(n s ). 

Therefore, from each considered partition we obtain the contribution 



<(vr)e n (7r) 

lim — = uJ-K)d v , n , . . . d Pl . . 

n^oo n q n s q Pi(1) Pl(a) 

where the division by n q comes from the normalization of the partial trace ip q {n). 
Collecting contributions associated with all n G AfC m (pi, qi, ■ ■ ■ ,p m , q m ), we obtain 

lim iP q (n)(S Pugi (n) . . . S Pm>qm (n)) = b *g( n ) 

n— too ' * 

ir£AfC'? n (pi,qi,...,p m ,qm) 

where b*(n) = b* q (n, f) corresponds to the unique coloring / of n defined by the tuple 
(pi, q±, . . . ,p m , q m ) and to the matrix B = DU . Finally, the proof that the expression 
on the right-hand side is equal to tp q (u Pltqi . . . ojp myqm ) is similar to that of Lemma 5.1, 
which completes the proof. ■ 

Corollary 7.1. Under the assumptions of Theorem 7.1, the joint ip{n)- distributions of 
blocks and block units converge to the joint tp- distribution of the truncated matricially 
free Gaussian operators and truncated units, respectively, as n — > oo. 

Proof. If we replace ip q (n) by vp(n) at the end of the proof of Theorem 7.1, we need 
to sum over q G [r] the corresponding right-hand sides with weights d q , respectively, 
which result from normalizations of partial traces. Therefore, in order to obtain a 
combinatorial expression for the corresponding mixed moment, it suffices to replace 
b* q {j) in the above formula by b*(n) = J2 q dqb*(ft)- This gives ip(uu Puqi . . . 0J Pmtqm ), which 
proves our assertion. ■ 

Corollary 7.2. Under the assumptions of Theorem 7.1, the array (S Pjq (n)) is asymp- 
totically matricially free with respect to the array of states (tp Piq (n)) on A(n) defined by 
the family of normalized partial traces {ip q {n))i< q < r as n — >■ oo. 

Proof. It follows from Theorem 7.1 that the T/> 9 (n)-distribution of (S Ptq (n)) converges 
to the ^-distribution of (ou Ptq ) for any p, q. This, Theorem 7.1 and Proposition 4.2 give 
our assertion. ■ 

Theorem 7.2. Under assumptions (A1)-(A4), the joint (p{n)- distributions of blocks 
and block units converge to the joint ip- distribution of the matricially free Gaussian 
operators and the associated units, respectively, as n — > oo. 

Proof. The proof is similar to that of Theorem 7.1 and is based on the combinatorial 
arguments of type used in [5, Lemma 6.1] which lead to [5, Lemma 6.2]. ■ 

Corollary 7.3. Under assumptions (Al)-(AJf.), the array (S Pjq (n)) is asymptotically 
matricially free with respect to the array ((p p , q (n)) of states on A{n) defined by <p(n) 
and the family of normalized partial traces (ip q (n))i< q < r as n — >■ oo. 



Proof. The assertion follows from Theorem 7.2 and Proposition 4.1. 
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8. Symmetric matricial freeness 

In order to compare the asymptotics of blocks of random pseudomatrices with that 
of symmetric random blocks, we shall now introduce a symmetric analogue of matricial 
freeness. 

Roughly speaking, in this concept one replaces ordered pairs by two-element sets 
and assumes that the array (Aij) of subalgebras of a unital algebra A contains the 
diagonal and is symmetric. Moreover, each algebra Aij contains an internal unit l £ j 
which agrees with lj ti for any G J. By X we denote the unital algebra generated 
by the internal units and we assume that it is commutative. By ((pij) we denote an 
array of states on A. 

Instead of sets A m , we shall use their symmetric counterparts, namely subsets of I m 
of the form 

n m = {{{h, i 2 }, {i 2 , h}, {im, im+i}) ■ {h, h) ^ {h, k} ^ ■ ■ ■ ^ {im, Wi}} 

where m G N, with their union denoted II = Um=i ^m- ^he mam difference between 
'symmetric matricial freeness' and matricial freeness is that in all definitions we have 
to use II instead of A. 

Definition 8.1. We say that the array (ljj) is a symmetrically matricially free array 
of units associated with (Aij) and ((pij) if for any diagonal state ip it holds that 

(1) ip(uiau 2 ) = <f(ui)ip(a)ip(u 2 ) for any a G A and Ui,u 2 G X, 

(2) if a k G A ik j k n Kerip ik j k , where 1 < k < m, then 

n , _ f <p(aa 2 ...a m ) if ({h, . . . , {i m , j m }) G II 

^«W«2..-a m )-| Q otherwise 

where a G A is arbitrary and {ii,ji} ^ ... ^ {i m ,jm}- 

Definition 8.2. We say that a symmetric array (Aij) is symmetrically matricially free 
with respect to ((pij) if 

(1) for any a k G Ker(p ikJk nA ikJk , where k G [m] and j ± } ^ . . . ^ {i m , j m }, and 
for any diagonal state tp it holds that 

if(a x a 2 ■ ■ .a m ) = 

(2) (ljj) is a symmetrically matricially free array of units associated with (Aij) and 
(Vi,j)- 

The array of variables (a^j) in a unital algebra A will be called symmetrically matri- 
cially free with respect to (<Pij) if there exists a symmetrically matricially free array of 
units (hj) m A such that the array of algebras (Aij), each generated by a^j + aj ti and 
lij, respectively, is symmetrically matricially free with respect to (<fiij)- The definition 
of *-matricially free arrays of variables is similar to that of *-matricially free arrays. 

We shall need the symmetrized operators on the matricially free-boolean Fock space 
T of Section 4, like the cD^-'s defined as 
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In a similar way we define Qj, Qj, pij, etc. All these arrays are associated with matrix 
A = (cxij), which is supressed in the notation. In turn, the symmetrized units on T 
and J 7 Cf2, are of the form 

lj,jr' = + — lijljji an( f = 

respectively, for any i,j, where (lij) is the array of canonical units on T and P stands 
for the projection onto J6Cfi. More explicitly, is the projection onto the subspace 
of T spanned by tensors which begin with or e^k for some k, and, in addition, by 
Q if % = j. 

Proposition 8.1. If the matrix A = (etij) is symmetric, then 

(1) the ipj- distribution of Qij is the semicircle law of radius 2ai,j for any 

(2) the array (fiij) is symmetrically matricially free with respect to (ipij)- 

Proof. lfi = j, then the first claim is similar to (1) of Proposition 4.1 since the action 
of Ujj onto ejj is similar to the action of Qj onto the vacuum vector. If i ^ j, then for 
even and positive m, it holds that 

^K}) = E E ^ (^iun) ■ ■ -t^jj) 

ei...,e m 6{l,*} (h,ji),-,(i m ,jm)e{(i,j),(j,i)} 

and there is a bijection between NC 2 m and products of the form ^(e^jj . . . i em (e im j m ) 
whose action onto ejj is non-trivial. This bijection is obtained as follows. If n e NC 2 rn 
is given and {l,r} is a block of n, where / < r, then e; = *, e r = 1 and = (i r ,j r ) 

with 

/. ■ \ _ f (jo(r),io(r)) if has an outer block 
tV ' Jrj ~ \ ' ' otherwise 

where o(r) is the right leg of the nearest outer block of {l,r}. It can be seen that 
this mapping is onto since in order to get a non-trivial action on ejj of the considered 
product of operators, the action of each £(eij) (corresponding to the right leg of some 
block) must be followed by the action of l{e^i) (corresponding to the right leg of another 
block) or £*{eij) (corresponding to the left leg of the same block). In particular, acting 
first with £(eij) on a vector from T and then with £{e^j) or £*(ejj) gives zero. This 
completes the proof of (1) since the even moments of the semicircle law of radius 2aij 
are given by M m = a™c m / 2 , where c m / 2 = |A/"C^|, m G 2N, are Catalan numbers. The 
proof of (2) is similar to that of (2) of Proposition 4.1. We show the slightly more 
general result that the array Bij = C(p i j -, p*j, t^) is symmetrically matricially free 
with respect to (ipij), where we use relations 

"""""" _ /-. -I- 

and normalization ipj(tij) = 1 for any (i,j) G J. The details are left to the reader. ■ 

For a given partition [n] = N 1 U N 2 U . . . U N r of the form considered before, it is 
useful to introduce the notation 

N p , q = (N p x N q ) U (N q x N p ) 
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for sets of pairs which label 'symmetric blocks' of random pseudomatrices. If we let 
n — > oo and assume that blocks grow as before, we can obtain the asymptotic behavior 
of these symmetric blocks. 

Theorem 8.1. Under the assumptions of Theorem 7.1, the array of symmetric blocks 
given by 

is asymptotically symmetrically matricially free with respect to (ip p<q (n)) as n — > oo. 

Proof. Observing that 

z ( n ) = ! S P,p( n "> X P = q 

p ' 9 [ Sp,q( n ) + Sq,p( n ) otherwise ' 

it suffices to use Theorem 7.1 and Proposition 8.1 to prove the assertion. ■ 

9. Symmetric random blocks 

We show in this section that the tracial asymptotics of random pseudomatrices is the 
same as the asymptotics of complex Gaussian random matrices and that this similarity 
can be carried over to the level of their symmetric random blocks. 

The context for the study of random matrices originated by Voiculescu [10] is the 
following. Let \x be a probability measure on some measurable space without atoms and 
let L = ni< p <oo ^(/-O be endowed with the state expectation E given by integration 
with respect to fi. The *-algebra of n x n random matrices is M n (L) = L<g)M n (C) with 
the state r(n) = E <g> tr(n), where tr(n) is the normalized trace. 

In order to compare the asymptotics of the blocks of random matrices with that 
of random pseudomatrices, we partition each set [n] into disjoint non-empty intervals 
as in the case of pseudomatrices and we set again D = diag(<ii, d 2 , . . . , d r ) to be the 
associated diagonal dimension matrix. 

By a complex Gaussian random matrix we understand a matrix (Yi,j(n))i<ij< n , in 
which Yij{n) = Yj :i (n) for any i,j,n and 

{ReF M (n)|l < i < j < n} U {ImF M (n)|l < i < j < n) 
is an independent set of Gaussian random variables. Its submatrices of the form 

T P,q( n )= Yl Y i:j (n) ® eiJ (n) 

will be called symmetric random blocks, where p, q G [r], and {ejj(n)|l < i < j < n} is 
a system of matrix units. 

We will study the asymptotics of symmetric random blocks under some natural as- 
sumptions. Our setting is very similar to that in [10,12] except that we will assume 
that the variances of |Yy(n)|, where i,j G [n] and n is fixed, are block-identical rather 
than identical. The case of non-Gaussian random matrices [2] and the corresponding 
symmetric blocks can be treated in a similar way. 

In order to find a Hilbert space realization of the limit joint distribution of symmetric 
random blocks under r(n), we will use the symmetrized truncated Gaussian operators 
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u p>q on T and the state ip — J2j djipj on B^F), where the dependence of tu Ptq s on the 
matrix B = DU is supressed in the notation. 

Theorem 9.1. Let (^ / i,j(^))i<ij<n be a complex Gaussian random matrix for eachn G N 
such that 

(!) E ( Y i,j( n )) = f or an V h3i n , 

(2) E(|Yjj(n)| 2 ) = u Ptq /n for any i G N p and j G N q and any n, 
where U := (u Piq ) G M r (M). Then it holds that 

lim T(n)(T Puqi (n) . . .T Pmtqm (n)) = 4>(Q Pl , qi . . .u) Pm , qm ) 

n—toc 

for any pi, ?i, . . . ,p m , <? m G [r], where uj p>q 's are associated with B := DU, where D is 
the dimension matrix. 

Proof. Our proof refers to the original proof of Voiculescu [10,12], which is followed by 
some combinatorial arguments referring to non-crossing pair partitions (in our approach 
the variances may vary). Let r q (n) = E Cg> tr g , where tr q (A) = l/n q ^2j eN<i Ajj is the 
normalized partial trace of A G M„(C) corresponding to q G [r]. We have 

T q (n)(T puqi (n)...T Pm , qm (n)) 
Yl E ( Y h,h( n ) ■ ■ ■ Y imjM) tr q (e iujl . . . e imtj J. 

{h,ji)eN P1}qi ,...,(i m ,j m )eN Pm}qm 

Since all terms are zero for m odd, throughout the rest of the proof we assume that 
m = 2s for some s G N. 

An individual term of this sum is then non-zero only if 

ji = h,h = h, ■ ■ ■ , 3m = h e N q 

and there exists a bijection 7 : [to] — > [m] such that 

7 2 = id, 7(fc) ^ k and i 7(fc) = j 7(fe) = i fc . 

Each such term is 0(n _s_1 ) as n — > 00. For convenience, we restrict our attention to 
the case when 

(ii, ji) G iV Pl xJV, ir .., (i m , j m ) G A^ Pm x AT 9m , 

since the remaining cases can be treated in a similar way, with some pi's interchanged 
with the corresponding q^s. 

The number of non-zero terms of the considered type is ©n(7), where 7 runs over 
the set of permutations of [to] such that 

7 2 = id, 7(fc) 7^ fc and p fe = g 7(fc) , g fc = p l{k) 

for any G [to]. If q — pi, we obtain 

©n(7) = card{(ii, . . . , i m ) G N pi x . . . x N Pm : i fe = i^k)+i, ^fc+i = ^ 7 (fc)}, 

where addition is modulo to. Denote by 71(7) the pair-partition of [to] defned by 7. 

If 77(7) is a crossing pair-partition, then 0„(7) < 0(n s ), as in the free case, which 
will make the corresponding mixed moments disappear as n — > 00 since we still need 
to multiply these cardinalities by the corresponding expectations, which are 0(l/n s ), 
and divide them by n q due to the normalization of the partial trace. 
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<2) 



i(4) 



l(2r) 



inner blocks inner blocks ... inner blocks 

m(l)n(l) n(2) ra(3) n(4) n(2r-l) n(2r) m(2) 



Figure 2. Inductive step involving subpartitions of 77(7). 



In turn, if ^(7) G A/"C^ \ A/"C^({pi, qi}, {p m , q m }) or q ^ p 1 , then again n (7) < 
0(n s ). On the other hand, if 71(7) e jVC^({pi, 51}, ... , {p m , g m }) and g = pi, then the 
number of independent indices which enter in the computation of n (7) is s + 1. We 
claim that for such 7 we have 



e„(7) = card{(ii,i r( i), . . . ,i rW ) G N Pl x N Pr(1) x . . . x A^ Pr(s) 



where 



^(7r( 7 )) = {r(l),r(2),...,r( S )}. 

In fact, it is not difficult to show that the conditions which define O n (7) can be reduced 
to s + 1 independent indices ii,i r (i), ■ ■ ■ ,i r (s) which can be interpreted as independent 
colors assuming arbitrary values from the corresponding intervals N Pl , N Pr(1) , . . . , N p ,. , 
respectively, with index %\ coloring the imaginary block. 

The proof of this fact essentially reduces to the inductive step involving subpartitions 
of 77(7) of the form shown in Fig. 2, where the covering block can also be interpreted as 
the imaginary block associated with 77(7). We want to show that the conditions which 
define 8,1(7) reduce all indices involved here to those associated with the right legs of 
the blocks of the subpartition shown above. In particular, on the level of blocks of 
depth one, we have pairings 

7(n(l)) = n(2), 7 (n(3)) = n(4), . . . , 1 (n(2r - 1)) = n(2r) 

which, in view of the condition i k = i 7 ( fe ) +1 , lead to the equation 



(3) = 



= ^ J 



(2r-l) — 1m{2) 



which says that the indices associated with the left legs of blocks of depth one are equal 
to the index associated with the right leg of its covering block. No blocks of depth 
greater than one lead to new conditions involving i m ( 2 ) and for that reason the index 
i m ( 2 ) can be used to color the covering block of the considered subpartition. The same 
pattern is repeated for blocks of arbitrary depth which are covered by the same block. 
The same holds for blocks of zero depth, where the role of the covering block in the 
above reasoning is played by the imaginary block (in that case the imaginary block is 
colored by 

This argument can be viewed as a proof of the inductive step of our claim that the 
conditions involving the indices i±, i 2 , ■ ■ ■ , i m , used to define the numbers 6 n (7) for the 
'right' 7 reduce to s + 1 independent indices associated with the right legs of the blocks 
of 71. If the covering block is interpreted as the imaginary block, we obtain the starting 
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case of the induction. This proves our formula for O n (7), from which we easily get 

©nw) = n Pi n Pr(l) ■ ■ ■ U Pr(s)- 

Let us remark here that not all tuples i r (i), ■ ■ ■ , i r (s)) which contribute to O n (7), 
but only those which are pairwise different, are actually used in our computation of 
the asymptotic joint distribution. If any two (or more) of these indices coincide, the 
corresponding contribution from the associated expectations is 0{n~ l ) and therefore 
becomes irrelevant in the limit. 

For those independent indices which are pairwise different and correspond to parti- 
tions 71(7) G J\fC 2 m ({pi, gi}, . . . , {Prn,Qm}), the corresponding expectation of Gaussian 
random variables K(Y il j 1 (n) . . . Y im j m {n)) is the product of 

E (^(™)^ 7(fc) j 7W (™)) = v ikJk (n) 

where j(k) G 7^.(77(7)), with i 7 ( fe ) = jk and j^) = ik- Now, since ik = i (k), the above 
expectation is equal to Vij(n), where % = i (u\ and j = In other words, j is the color 
assigned to block {k,j(k)} since we have shown before that we color the blocks with 
indices which correspond to their right legs and i is the color assigned to its covering 
block. This includes the case of blocks which do not have an outer block, then the role 
of the latter is played by the conditional block colored by ii = j m . 

Now, using our assumption on the variance matrices, Vij(n) = Vj^(n) = u qjP /n when- 
ever i G N p , j G N q . Therefore, taking into account all pairings and using the definition 
of symbol Uj(n, f) of Section 2, we obtain 

m i M...Y imJ Jn)) = ^fl 

for the relevant Gaussian expectations of variables associated with 71(7), where / G 
F r (7r(7)) is the coloring of the associated 71(7) defined by indices p r (i),p r (2), ■ ■ ■ ,p r (s) 
assigned to the blocks of 71(7), with color p 1 assigned to the imaginary block. 

Collecting expectations corresponding to all ix G MC 2 m and taking into account that 



Bm «^ = ^ (i| ...^ () , 



we obtain 



lim r q (n)(T puqi (n) . . .T PmAm {n)) = 

ireAfC^ ({pi ,qi},...{p m ,9m}) 

where b q (7r) = b q (-K, f) is the sum over all admissible colorings of ir and each b q (7r, f) 
coresponds to the matrix B = DU and to the color q of the imaginary block. 

If we replace T q (n) by r(n) in the above expression, we just need to sum over q G [r] 
the corresponding right-hand sides with weights d q , respectively, which result from 
normalizations of partial traces. Therefore, in order to obtain the required expression 
for the corresponding mixed moment, we need to replace b q (n) by b(ii) = J2 q d q b q (n), 
which completes the proof of the combinatorial formula. 

It remains to show that the mixed moment of the symmetrized Gaussian operators 
in the state tp gives the apropriate sum of combinatorial expressions obtained above. 
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This, in turn, is similar to the proof of Lemma 5.1. Therefore, our proof is completed. ■ 

Remark 9.1. In this context, let us remark that the array of symmetric random blocks 
(T pq (n)) is not symmetrically matricially free with respect to (r P)9 (n)) for finite n. 
In order to see this, it suffices to compute some simple examples of mixed moments 
of the T Pi(3 (n)'s associated with crossing partitions. For instance, if n — 3 and the 
blocks are one-dimensional, then T^iTi^T^T^Tx^T^T^i) ^ 0, where T\ = ti(3) and 
T PtQ = T Ptq (3) for any 1 < p, q < 3, which shows that condition (1) of Definition 8.2 is 
not satisfied. 

However, in view of Proposition 8.1 and Theorem 9.1, we can expect that symmetric 
random blocks are symmetrically matricially free as n — > oo. The notion of asymptotic 
symmetric matricial freeness is analogous to that of asymptotic matricial freeness (see 
Definition 7.1). It remains to define appropriate arrays of block units. Thus, let 

j€N p UN q 

for any p,q G [r] and each n G N. It is easy to see that is an internal unit in the 

algebra generated by T Pjq (n) and l Pjq (n) for any given p, q G [r] and n G N. 

Theorem 9.2. Under the assumptions of Theorem 9.1, let (j Ptq {n)) be the r ■-dimensio- 
nal array of states on the algebra A(n) = L ® M n (C) defined by the normalized partial 
traces r q (n), where q G [r], for any n G N. Then the array of symmetric random blocks 
(T Pt q(n)) is asymptotically symmetrically matricially free with respect to {r Ptq {n)). 

Proof. In view of Proposition 8.1, it suffices to show that the r m (n)-distributions 
of blocks T p>q (n) and units l Ptq (n) converge to the ?/> m -distribution of the symmetrized 
Gaussians Q Piq and summetrized units as n — > oo for any given m. It suffices to 
verify that l Ptq (n) leaves invariant any vector of the form 

wi(T Pli9l (n)) . . .w k (T Pktqk (n))e m 

where the product of polynomials Wi, . . . , Wk is in the symmetrically matricially free 
kernel form with respect to the array (r P)9 (n)) and {p, q}r\{p±, q±} ^ and kills all other 
vectors. This is quite easy to observe since any vector of this type is a linear combination 
of the form ^ i&Np cxiai <g> C{ + ^jgAr 9 Pjbj <S> e-,-, where a^, bj G L and ctj, j3j G C for any 
and thus it is left invariant by l Pjq (n), where (e,) is the canonical basis in C n . ■ 

10. Asymptotic freeness and asymptotic monotone independence 

If we make additional assumptions on the variance matrices (v it j(n)) of matricially 
free arrays (JQj(n)), we obtain asymptotic freeness which refers to the rows of random 
pseudomatrices. The proposition given below can also be viewed as an operatorial 
version of [5, Proposition 8.2]. 
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Proposition 10.1. If the arrays (Xjj(n)) of Theorem 3.2 are square and have identical 
variances within blocks and rows, then the sums 

r 

S p (n) := S Pj g(n), where 1 < p < r, 

9=1 

are asymptotically free with respect to both tp(n) and ib(n) as n — >■ oo. 

Proof. In view of Theorems 7.1-7.2, it suffices to show that the variables £i, . . . , ( r 
are free with respect to ip, and that oui, . . . , ui r are free with respect to ip, where 

r r 

(p = Y1 ^ and u p = Y1 

q=l q=l 

with the notations of Section 4. If {ei, e2, . . . , e r } is an orthonormal set of vectors, we 
have the natural isomorphism 

re: jr(0Ce g ) 

9 

of Remark 4.1 and then, since Q"p^q — 

for any p, q G [r] , we have 

where u>(e p ) = £(e p ) + f (e p ) is the canonical free Gaussian operator on J^i^&qCeq) for 
any p. Let us remark that u p and u>(e p ) denote different operators in our notation. 
Since w(ei), . . . , u(e r ) are free with respect to the vacuum state on B(J z ((§) q Ce g )), the 
operators Ci, • • • , ( r are free with respect to the vacuum state on B^F). Similarly, 

w p = a PiP 7*w(e p )7 

for any p G [r], where 7 = Pok and P is the canonical projection from J-(@ q Ce q ) onto 
the orthocomplement of CQ. Now, since Pu(e\)P, . . . , Pcu(e r )P are free with respect 
to VvC) = (- e qi e q) f° r an y 9) the operators Wi, . . . ,u r are free with respect to tp q for 
any q. Moreover, the ?/>p-distribution of u p does not depend on p since the variances are 
assumed to be identical in each row. Therefore, u>i, . . . ,u r are free with respect to ip. 
This completes the proof. ■ 

Analogous results hold for block-triangular arrays and lead to monotone indepen- 
dence (lower-block-triangular arrays) and anti-monotone independence (upper-block- 
triangular arrays). One has to remember that order is important for these notions of 
independence and therefore in that case we use sequences of variables instead of fam- 
ilies. We formulate the result only for lower-triangular arrays since the case of upper- 
triangular arrays is completely analogous. This result can be viewed as an operatorial 
version of [5, Proposition 8.3]. 

Proposition 10.2. If the arrays (X it j(n)) of Theorem 3.2 are lower-block-triangular 
and have identical variances within blocks and rows, then the sums 

v 

S p (n) := S Pt g(n), where 1 < p < r, 
9=1 
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are asymptotically monotone independent with respect to (pin) as n — >■ oo. 

Proof. Let ( p = Y7 q =i Cp,g> where p G [r]. In view of Theorem 7.1, it suffices to show 
that the operators d, . . . ,(r are monotone independent with respect to (p. Let T\ be 
the subspace of ^"(0^ Ce g ) of the form 

r 

m=l pi>...>p m fii,...,n m eN 

where H g = Ce g for any q G [r] and {ei,e2, . . . ,e r } is an orthonormal basis in some 
Hilbert space and denote by Q : -7 7 (® (J Ce q ) — > T\ the corresponding canonical projec- 
tion. Observe that we have 

C P = a PtP (3*ui(e p )(3 

where = Q o k and k is the same as in the proof of Proposition 10.1. Since the op- 
erators Qu(ei)Q, . . . , Qu(e r )Q are monotone independent with respect to the vacuum 
state on the free Fock space, the operators £i> . . . ,(r are monotone independent with 
respect to (p. This completes the proof. ■ 

Finally, we obtain asymptotic boolean independence of blocks of block-diagonal pseu- 
domatrices. 

Proposition 10.3. If the arrays (Xij(n)) of Theorem 3.2 are block- diagonal and have 
identical variances within blocks, then the sums S P)P {n), where 1 < p < r, are asymp- 
totically boolean independent with respect to (p{n) as n — )■ oo. 

Proof. This is a simple consequence of Theorem 7.2 and the fact that the family 
i(p,p)i<p<r is boolean independent with respect to the vacuum state <p on J 7 . ■ 
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